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Abstract 

We study the domain walls connecting different chirally asymmetric vacua 
in supersymmetric QCD. We show that BPS - saturated solutions exist only 
in the limited range of mass m < m* ~ 0.8| < Tr When m > m*, 

the domain wall either ceases to be BPS ~ saturated or disappears altogether. 
In any case, the properties of the system are qualitatively changed. 


1 Introduction 


Supersymmetric QCD is the theory involving a gauge vector supermultiplet V and 
a couple of chiral matter supermultiplets belonging to the fundamental repre¬ 
sentation of the gauge group SU{N). The lagrangian of the model reads 

C = ^Tr Jd^ew^ + d^ed^e J 

where Sf = efgS^ (for further conveniences, we have changed a sign of mass here 
compared to the standard convention). The dynamics of this model is in many 
respects similar to the dynamics of the standard (non-supersymmetric) QCD and, 
on the other hand, supersymmetry allows here to obtain a lot of exact results ]l|. 

Like in the standard QCD, the axial Da( 1) symmetry corresponding to the chiral 
rotation of the gluino held and present in the tree-level lagrangian (|1.1D is broken 
by anomaly down to Z 2 n- This discrete chiral symmetry can be further broken 
spontaneously down to Z 2 so that the chiral condensate < Tr > is formed. There 
are N different vacua with different phases of the condensate 

< Tr fc = 0,..., iV - 1 (1.2) 


It was noted recently that on top of N chirally asymmetric vacua (|1.2|) , also a 
chirally symmetric vacuum with zero value of the condensate exists. 

The presence of different degenerate physical vacua in the theory implies the 
existence of domain walls — static held conhgurations depending only on one spatial 
coordinate {z) which interpolate between one of the vacua at z = — 00 and another 
one at z = cx) and minimizing the energy functional. As was shown in ^j, in many 
cases the energy density of these walls can be found exactly due to the fact that 
the walls present the BPS-saturated states. The key ingredient here is the central 
extention of the N = 1 superalgebra P, ^ 


Wi«U=4(3).,W Av 






A domain wall is a conhguration where the integral of the full derivative in the RHS 
of Eq. (|1.3|) is non-zero so that the standard N = 1 SUSY algebra in the wall sector 
is modihed. A BPS-saturated wall is a conhguration preserving 1/2 of the original 
supersymmetry i.e. a conhguration annihilated by the action of two certain real 
linear combinations of the original complex supercharges Qa []• 

Combining (pT^) with the standard SUSY commutator {Qa,Q^} = 2((J^)fP^ 
and bearing in mind that the vacuum expectation value of the expression in the 

^Such BPS-saturated walls were known earlier in 2-dimensional supersymmetric theories (they 
are just solitons there) and were considered also in 4-dimensional theories in stringy context 

§• 
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square brackets in Eq. (|1.3| ) is zero due to Konishi anomaly [0 0, it is not difficult to 
show that the energy density of a BPS-saturated wall in SQCD satisfies a relation 


e = 


N 

Svr^ 


< Tr >oo - < Tr 


(1.4) 


where the subscript ±cxo marks the values of the gluino condensate at spatial infini¬ 
ties. 


The relation (|1.4|) is valid assuming that the wall is BPS-saturated. However, 
whether such a BPS-saturated domain wall exists or not is a non-trivial dynamic 
question which can be answered only in a specific study of a particular theory in 
interest. 

In Ref. [0 this question was studied for the SU{2) gauge group in the framework 
of the effective low energy lagrangian due to Taylor, Veneziano, and Yankielowicz [§. 
The situation is particularly simple when the mass parameter m in the lagrangian 
( |lI1) is small compared to Asqcd = A. In this case, chirally asymmetric vacua 
are characterized by large values of the matter scalar field x- The theory involves 
two different energy scales, and one can integrate out heavy fields and to write the 
Wilsonean effective lagrangian describing only light degrees of freedom. It is the 
lagrangian of the Wess-Zumino model with a single chiral superfield X and the 
superpotential 


3X2 


^X^ 

2 


(1.5) 


The corresponding potential U = \dW/dx\^ has two different non-trivial minima 
at < = ±^j4:A^/3m. A domain wall interpolating between these 

vacua is BPS-saturated. The solution can be found analytically [Q 


X(z) = X, 


X jQ‘irn,{z-zo) 
y/l g8m(z-zo) 


( 1 . 6 ) 


where Zq is the position of the wall center. 

In this approach, we are not able, however, to detect and study a chirally symmet¬ 
ric vacuum with < x >= 0 and the corresponding domain wall. Chirally symmetric 
vacuum appears when taking into account also the degrees of freedom associated 
with the gluon and gluino helds. The full TVY effective lagrangian is, again, a 
Wess-Zumino model involving two chiral superfields $ and X with the superpoten¬ 
tial 


W 



In 


$3x2 



(1.7) 


The corresponding potential for the lowest components 0, x of the superhelds <h, X 


m-p.x) = 


dW 

2 

dW 

2 


2 

( V\ 

d^p 

+ 

dx 

= 4 


+ 

3xV 


( 1 , 8 ) 


^See recent g for detailed pedagogical explanations. 
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(in the remainder of this and in the next section we will measure everything in the 
units of A) has three non-trivial minima: 


= A = 0 


(1.9) 


= e 


/3m\ 

t/3 


1/4 

\3mJ 


(±\ 


X = 


'(-) 

\3mJ 


1/4 


( 1 . 10 ) 


There are also the minima with inversed sign of (j) and y, but they are physically the 
same as the minima (|1.10|) : the vacuum expectation values of the gauge invariant 
operators < Tr > = (327r^/3) < > and < s^Sf > = < > (s-f is the squark 

held) are the same. 

The effective theory ( |1.5|) is obtained when the heavy degree of freedom 0 is 
frozen in the Born-Oppenheimer spirit so that = 1- the opposite limit 

m 1, we can freeze instead the heavy matter helds = 40^/3m and arrive 
at the Veneziano-Yankielowicz effective lagrangian [|T^ for pure supersymmetric 
gluodynamics which involves only the held 0.) 

Generally, one should study the theory with the potential (|1.8| ). The status of this 
ehective theory is somewhat more uncertain than that of (|1.5|) — for general value 
of mass, the TVY ehective lagrangian is not Wilsonean; light and heavy degrees of 
freedom are not nicely separated. But it possesses all the relevant symmetries of 
the original theory and satishes the anomalous Ward identities for correlators at 
zero momenta. We think that the use of the TVY lagrangian is justihed as far as 
the vacuum structure of the theory is concerned. 

In Ref. [Q the domain walls interpolating between the chirally symmetric min¬ 
imum (HU) and a chirally asymmetric minimum in Eq.( l.lU|) were studied along 
these lines. It was shown that these walls are BPS-saturated at any value of mass. 

In this paper, we study the complex domain walls interpolating between diherent 
minima in (|1.10|) . Rather surprisingly, we have found that the solution of the BPS 
equation exists only for small enough masses m < m* = 4.6705 .... At larger values 
of mass, the BPS equations have no solution. A phase transition occurs. 


2 Solving BPS equations. 

BPS equations for the domain wall in a generalized Wess-Zumino model with two 
chiral superhelds read P, ||, ^ 

dz(p = ±dW/d(j), dzX = i^dW/dx (2.1) 

^To see that, one should use the amended potential of Ref. Q which is “glued” of different sectors 
related to the different branches of the logarithm; cf. an analogous situation in the Schwinger model 
|pd|. For our present purposes, these complications are irrelevant, however. 
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In our case, the superpotential is given by the expression (pTTI). Let us choose the 
positive sign in Eg . (12.11) and try to solve it with the boundary conditions 


(j){-oo) = , (j){oo) = e~ 

/ 4 \ l/4 


t/3 (^\ 


V 4 y 


4 


x(-oo) = f-^') , x(oo)=^fA) 


\3mJ 


\3mJ 


( 2 . 2 ) 


(the negative sign in ( [2.1D would describe the wall going in the opposite direction in 
z). The solution of the equations ( p.l|) with the boundary conditions (|2.2|) has the 
hxed energy which coincides with (|1.4|) after the proper normalization 


327r2 


Tr 


(2.3) 


is chosen. 

Technically, it is convenient to introduce the polar variables y = 0 = 

Then the system (p.l|) (with the positive sign chosen) can be written in the form 


dzP = —mpcos(2a) + ^ cos(3/3) 
dzOt = msin(2a;) — sin(3/9) 

dzR = 1R? [cos(3/d) — sin(3/3)(3/d + 2a)] 

dzl3 = —2i? [sin(3/9) ln(i?^p^) + cos(3/d)(3/9 + 2a)] 


(2.4) 


The wall solution should be symmetric with respect to its center. Let us seek for 
the solution centered at 2 ; = 0 so that 

p{z) = p{—z)^ R{z) = R{—z), a{z) = 7r/2 — a(— 2 ;), (3{z) = —7r/3 — (3{—z) (2.5) 
Indeed, one can be easily convinced that the Ansatz ( p.5| ) goes through the equations 

(ID- 

The system (p.l|) has one integral of motion [^: 


Im W{(f), x) = const 


( 2 . 6 ) 


Indeed, we have 




d(j) 


dx 


dW 


d(f> 


+ 


dW 


dx 


= dM 


It is convenient to solve the equations ( p.4|) numerically on the half-interval from 
2 ; = 0 to 2 ; = cx). The symmetry (p^.5|) dictates a(0) = 7r/4, /9(0) = —n/6. The 
condition (|2.6|) [in our case Im W (0, x) = 0 due to the boundary conditions (|2.2|) ] 
implies 




\n{R^p‘^) — 1 + mp^ = 0 


(2.7) 
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Figure 1: Mismatch parameter A as a function of mass 

Thus, only one parameter at z = 0 [say, i?(0)] is left free. We should ht it so that 
the solution would approach the complex minimum in Eq. (|1.10[ ) at z ^ oo. 

It turns out that the solution of this problem exists, but only in the limited 
range of m. If m > m* = 4.6705..., the solution misses the minimum no matter 
what the value of 7?(0) is chosen. This is illustrated in Fig. 1 where the “mismatch 
parameter” 

A = min min \l\x{.z) — + \(l){z) — (2.8) 

R{0) ^ 

is plotted (in a double logarithmic scale) as a function of mass. The dependence 
A(m) hts nicely the law 


A(m) = 0.56(m — m*) 


0.44 


(2.9) 


It smells like a critical behavior but, as A is not a physical quantity, we would not 
elaborate this point further right now. 

At m = 4.6705 or at smaller values of mass, the solution exists, however. The 
prohles of the functions p{z) and R{z) for three values of mass: m = 0.2, m = 2.0, 
and m = 4.6705 in the whole interval —oo < z < oo are plotted in Figs. 2,3. For 


small values of m, the solution approaches, as it should, the analytic solution ( p..6|) 
(with (j) = found in Ref.[Q, so that p{z) and R{z) stay constant. 

In Fig. 4 we plotted the dependence of R(0) on m. We see that 7?(0) has a 
singular behavior at the phase transition point, but stays non-zero. 
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Figure 4: The ratio rj = R{0)/R{oo) as a function of mass 


3 Discussion. 

We have found that the properties of the system are drastically changed at m = 
4.6705... A. It makes sense to express the result in invariant terms and to trade A 
for an invariant physical quantity such as the gluino condensate S = | < Tr > | 
in a chirally asymmetric vacuum. From Eqs.(|1.10|), (|2.3|), we obtain 


73 

Thus, the phase transition occurs at 

m* 


(3.1) 


(3,2) 


The particular numerical value ( |3.2|) was obtained by studying the effective TVY 
lagrangian with the potential (|1.8|) and the standard kinetic term + \dx\^- We 
cannot claim that the phase transition in the theory of interest ( 0 ) would occur 
at exactly the same value of mass. We believe, however, that the value of m* in the 
supersymmetric QCD should be close to (|^) . 

We hasten to comment that it is not a phase transition of habitual thermody¬ 
namic variety. In particular, the vacuum energy is zero both below and above the 
phase transition point — supersymmetry is never broken here. Hence Eyaci"^) = 0 
is not singular at m = m*. 

Some similarities may be observed with the 2-dimensional Sine-Gordon model 
where the number of the states in the spectrum depends on the coupling constant 
f3 so that the states appear or disappear at some critical values of [3 0. May be a 
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more close analogy can be drawn with the N = 2 supersymmetric Yang-Mills theory. 
The spectrum of the system depends on the Higgs expectation value u = < Tr >. 
A study of the exact solution of the model due to Seiberg and Witten []^ displays 
the existence of a “marginal stability curve” in the complex w-plane [T^. When 
crossing this curve, the spectrum pattern is qualitatively changed. 

To understand better the physics of the phase transition in supersymmetric QCD, 
one has to study in more details what happens in the region m > m*. There are 
two logical possibilities: 


• The complex domain wall solution may still exist at larger masses, but this 
solution cannot be BPS-saturated anymore. 0 In this case, it would be very 
interesting to study what happens in the limit m —>■ oo when the matter helds 
decouple and the theory is reduced to the pure N = 1 supersymmetric Yang- 
Mills theory. Do the domain walls interpolating between different chirally 
asymmetric vacua survive in this limit ? Q 

• The complex domain wall solution can disappear altogether at m > m*. 


The question of the existence and, if the solutions are there, of the properties 
of the complex domain walls at large masses is now under study. But in any case, 
it is clear now that there is no smooth transition between the weak coupling Higgs 
regime which is realized at small masses in chirally asymmetric phases and the strong 
coupling regime at large mass values. 
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